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Quark model matrix elements can be computed using bosonic operators and the holomorphic 
representation for the harmonic oscillator. The technique is illustrated for normal and exotic baryons 
for an arbitrary number of colors. The computations are much simpler than those using conventional 
quark model wavefunctions. 
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The non-relativistic quark model has been a useful tool 
in the study of hadrons. Baryons and mesons are de- 
scribed by quantum mechanical wavefunctions for non- 
relativistic constituent quarks. The lowest lying baryons, 
the 81/2 and 10 3 / 2 , are three quark states with wavefunc- 
tions which are completely antisymmetric in color, and 
completely symmetric in position and spin-flavor. 

The properties of baryons can be computed in a sys- 
tematic expansion in powers of 1 /N c , where N c — 3 is the 
number of colors in QCD. There is an exact 5C/(6) sym- 
metry for baryons in the N c —* oo limit of QCD, which 
is broken by 1/N C corrections pj. The systematic appli- 
cation of the 1/N C expansion requires computing pro 
erties of 5/7(6) baryon representations @>I3>0>La1 
S \Ei EH • The non-relativistic quark model also has 
SU(6) symmetry, and is a convenient way to keep track 
of SU (6) group theory. 

Applications of the quark model to compute baryon 
matrix elements, and to compute SU (6) Clebsch-Gordan 
coefficients, require evaluting quark operator matrix el- 
ements between baryon states. It is useful to have an 
efficient means of computing these quark model matrix 
elements. The traditional approach has been to write 
down color ® spin <g> flavor wavefunctions, and use them 
to compute matrix elements |l2t Il3j . The wavefunc- 
tion is completely antisymmetric in color, with the quark 
color indices contracted using e a bc- As a result, the three 
quarks must enter as one of each of the three colors, and 
the wavefunction is written by ordering the quarks by 
their color index. As an example, the proton (J z = 1/2) 
wavefunction is 

\p) = [uud\ ® [2 TU - TIT - ITT] + cyclic, (1) 

and the position of the quark in [uud] is the value of its 
color index. These wavefunctions are then used to com- 
pute operator matrix elements. For example, the matrix 
element 

(p\q^ 3 T 3 q\p) = (p\ -L f(l + 1 + 1) [uud\ ® [2 TTI] 



- (l - l - l) [uud] <g> [TIT] 

- (-1 + 1 - 1) [uud] [|TT] + cyclic 



where the terms in each parenthesis are the values of 
a 3 r 3 q for each of the three quarks. The above matrix 
element shows that in the quark model, the proton axial 
charge is (5/3)<jvt, where gA — 0.75 is the axial charge of a 
constituent quark 0] . Computations using quark model 
wavefunctions rapidly get complicated, particularly if one 
is interested in studying matrix elements for arbitrary N c , 
or for exotic baryons, both of which involve states with 
more than three particles. 

There is an alternate way to compute matrix elements 
which is more efficient, and more easily implemented for 
computer calculations. In the ground state baryons, and 
for exotic baryons, one is interested in hadron states 
which are completely symmetric in spin ® flavor. In this 
case, one can treat the quarks as bosons, completely ig- 
noring the color quantum numbers. One can then use 
bosonic operator methods developed for the harmonic os- 
cillator |15t flq to compute the matrix elements of color 
singlet quark operators, since they do not care about 
the color of individual quarks. The bosonic method has 
been usedpreviously to study the properties of large- N c 
baryons |J H, El > an d the relation between the quark 
model and chiral soliton model 0]- Some of the tech- 
niques used here are related to those in Ref. ^t|, where 
matrix elements in the N c — > oo limit were obtained us- 
ing spin- flavor coherent states. 

The one-dimensional harmonic oscillator can be quan- 
tized using creation and annihilation operators and a 
which satisfy the commutation relation [a, aH = 1. In 
the holo morp hic representation for the harmonic oscil- 
lator 0, Il6| . states are given as holomorphic functions 



ip(z) of a complex variable 
state 



z, a 



and a — d/dz. A 



|V.) - ^c n |n)=^^=(ot)«|0) 

n n Vn\ 



(3) 



in the occupation number basis corresponds to the func- 
tion 



i>{z) = £lM n = E-^T*" 



(4) 



10 ,. 5 

^ + cychc=-, 



(2) 



in the holomorphic basis, with expansion coefficients 
"0n = c n j \fn\. The inner product of two states is given 



2 



by 



00 f d 2 z 



-w 2 



[*(*)]* V(z), (5) 



where d 2 z/(27ri) = dRezdlmz. The operator d/dz is the 
hermitian conjugate of z using the inner product Eq. JSJ, 
as can be verified by integration by parts. 

The quark model bosonic algebra needed to study 
ground state and exotic baryons |19| is generated by 
quark and antiquark creation and annihilation operators 
q\ a , q la , q^ a , q IQ , where 1 = 1,2 is a spin index and a = 
1, . . . , F is a flavor index. In the holomorphic representa- 
tion, q| a and q^ la are multiplication by the complex vari- 
ables Q la and Q la , and q iQ ,q lQ are d/dQ la ,d/dQ la , re- 
spectively. Baryon wavefunctions are holomorphic func- 
tions ip (Qio, Q la ). Ground state baryons, which do not 
contain antiquarks, are functions ?/>(Q lQ ). An alterna- 
tive notation which is used also is Qn = Uf, Q21 = uj,, 
Q12 = d]i Q 11 = u, q 11 = Uf, etc. 1 Quark matrix el- 
ements can be computed using bosonic operators or the 
holomorphic representation, and the two methods are 
equivalent. The holomorphic representation is particu- 
larly convenient for computer implementation. 

Quark operators are differential operators in the holo- 
morphic representation. For example 



qWq 



d 
dut 



d 



dt 



d_ 

ddt 



d 



(6) 



Baryon wavefunctions are simpler than in the traditional 
approach. The wavefunctions of the J z — 1/2 octet and 
J z = 3/2 decuplet baryons with I = I3 are: 



(7) 





= 7I U T 


i u l d l 


- H d t] : 


, 1pA + + 


= 7e w T 
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- u i s t] . 
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- J_ s 3 



The matrix element of Eq. (j2J) is computed as 
qWq \p) = -^=u T [3w T d x + u x rf T ] , 
on applying Eq. JBJ to Eq. 0, so that 



(8) 



(pjqWqjp) = - [3 (u 2 d l {u^d^ ~ (u^u^ \u^ Ul d^ )] 



i [3-2-1] = -. 
3 L J 3 



The 9 + pentaquark is an exoticness E = 1 baryon, 
and is a qqqq IO1/2 state in the quark model. The wave- 
functions of the I = I3, J = 1/2 pentaquarks are 



2V3 
1 



x [2ut s 1«T + u l d l d i - UfdjJf - 2u x s T s T ] , 

10 

x [2tt T dj. J T + uj.s T s T - UT s i« T - 2ujd T J T ] , 



"3/2 



2V3 



(10) 



The wavefunctions of the 1=1 state 81 in the spin-1/2 
27, spin-3/2 27, and the I = 2 state 8 2 in the spin-3/2 
35 and spin-5/2 35 are 



^(27i /2 ,0i) 
^(27 3/2) 0i) 

y;(35 3/2 ,e 2 ) 
^(35 5/2 ,e 2 ) 



1 



2V3 
1 



u T [u T d x - ujd T ] [-?i T sx + uxs T ] 
U T 5 T I u T d l ^ u l d l] > 



2V2 
1 3 r 



30 

1 4- 

u T s T , 



2\/6 



(11) 



from which the other wavefunctions can be generated 
by applying raising and lowering operators. Matrix el- 
ements involving pentaquark states can be computed us- 
ing Eqs. man) . 

The qqqqq states in the quark model also include a 
81/2 • This qqqqq octet is a qqq octet plus a flavor singlet 
qq pair, and has exoticness E = even though it is a 
five-quark state. It is useful to have wavefunctions for 
baryons states containing qq pairs. Define 



S = qq = N- 



(12) 



which annihilates a qq pair, and its hermitian conjugate 
S' which creates a qq pair. We will need the commutation 
relations 

[5,5+] = 2N F + N q + N q , 
[A^.St] = A^-A^, 

[A^5] = -Af, 



[A^5] = At 



(13) 



(9) 



1 Q and Q are independent complex variables. The complex con- 
jugate of Q is Q* . Q is in the conjugate representation to Q, so 
that = Uf is a spin-up ft, and — u' = uy is a spin down S. 



where we use the notation of Ref. 19], = q* A A q, 
A A = qA A q, A A = q^A A q, where A A ,A A are spin- 
flavor matrices in the fundamental and anti-fundamental 
of SU(6). 

All minimal states satisfy the constraint that flavor 
singlet annihilation is impossible, and so obey 



51^ = 0. 



(14) 



3 



We now construct a non- minimal 8, denoted by 8', and 
defined by 



18') 



A 5+ 18) 



(15) 



where A is a normalization constant. The value of A is 
determined by computing the normalization: 



(8'|8') = |A| 2 <8 155^8) 

= |A| 2 (8|[5,5t]|8)=9|A| 2 , 

^A = i (16) 

In the holomorphic representation 

S^S = Q ta Q ia . (17) 

The matrix element of an operator such as the Hamil- 
tonian H can be computed similarly: 

(8'\H\T6) = A (8 \SH\ 10) = A (8 \[S,H]\ TO) , (18) 

and the commutator can be evaluated with the help of 
Eq. (| 1 31) . Matrix elements of AE = operators involving 
non-minimal states are equal to matrix elements of AE — 
— 1 operators involving the corresponding minimal state 
with the same flavor quantum numbers. Of particular 
interest are the matrix elements of T 8 , T 8 and T|, 

(8'|T 8 |T0) = 0, 
(8'|T 8 |To) = i(8|T 8 |T0), 

(8'|T 8 |T0) = -i(8|T 8 |T0), (19) 
o 

which are needed to analyze the masses of exotic baryons 
including SU(3) breaking [2tT |. 

The results so far have been for N c = 3. It is simple 
to extend the results to arbitrary N c . Let 



(20) 



which is a spin and isopspin singlet combination of u and 
d quarks, with 



(21) 



where Olq 1 "^ counts the number of u and d quarks. In 
the holomorphic representation 

The state 



(22) 



\b) = 



|0) , A = 



is a state with I = J 



I 3 = J 3 = 6/2. The state 



In, b) 



\b) , 1pn,b = 



\jb\o n ,b 



(23) 



(24) 



is a baryon with I = J = I 3 = J 3 = b/2 and N c = In + b 
quarks, and a n ^ is a normalization constant. One can 
similarly construct exotic baryons for arbitrary N c by 
multiplying the N c = 3 wavefunctions by powers of D n . 
The normalization a n 6 is given by 



<J n ,b = (b 



©«(S)t)' 



(25) 



This can be evaluated by expanding 2)™, (3D*)™ in a power 
series and summing the various terms. An alternate 
method is to use the commutation relation Eq. i|21[l . 
Commuting one 2) past (2)^) to annihilate \b) gives 



\b) 



r=0 
n-1 



]T (&Y + 2j (St 

r=0 

n(n + &+l)(2) t ) n_1 |6), 



n—r—l 



so that 



o" n ,6 = n(n + b+ l)cr n _i ib . 
Using do ,6 = 1 gives 

_ nl(n + b+ 1)1 
° n ~ (6 + 1)! • 



|6) 
(26) 

(27) 
(28) 



Eq. I|26l) can be derived more simply using the holo- 
morphic representation, and follows from the chain rule, 



1 j 1 ; du T dd L du L ddj V6! 



= n(n + b+l)D T 



(29) 



The matrix element of 2l 33 = qW 3 q can be computed 
similarly, using 



and 



2l 33 |6) = &|&), 

[2l 33 , S*] = 2) 1 , 
Mi + M T = 53, 



(30) 



2),2V 



= 21 



:s:i 



(31) 



Defining 



CTn,b€n,b — (b 

cr n ,br]n,b = < & 



2)™ 2l 33 (23 f ) 
23™ £ f (3)ty 



(32) 
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and commuting 2l 33 through in the first equation, and 
one 2) through in the second, as in Eq. Il2(ill gives the 
recursion relations 

= 2nr) nib + 6, (33) 
n(n + b + l)??„.b = (n- l)(n + b + 2)r) n -i,b + b. 

The same relations follow more simply in the holomor- 
phic representation using the differential form of the op- 
erator, Eq. JfJJ), and the chain rule. 

The solution to the recursion relations is 



Vn,b 
Cn.b 



6 + 2' 

b(2n + 6 + 2) _ b(N c + 2) 



(6 + 2) (6 + 2) ' 

using N c = 2n + 6. The matrix element for 6 = 1 is 

N r + 2 



9a(NNtt) = U,i = 



(34) 



(35) 



which is the value of the NNir axial coupling as a func- 
tion of N c in the quark model |2lJ, and reduces to 5/3 
for N c = 3. The matrix element of 2l 33 in the A as a 
function of N c is given by using 6 = 3, 



g A {AAir) = 



3(N C + 2) 



(36) 



One can similarly compute the N — > A transition axial 
matrix element of 2l ++ . The operator relations needed 
are: 

21++ = u\d L , [2t++,S) t ] = (u f ) 2 , 21++ |6) = 0, 



(ut)» = v /(6+l)(6 + 2)|6 + 2). 
By a method similar to the derivation of Eq. (|26I29(I 



(37) 



(n-l,6+2|2l ++ |n,6) 



1 /n(6+l)(n + 6 + 2) 



(6 + 3) 



1 / (&+l )(.V,-6l(-V,. + 6+ 4) 



(6 + 3) 



Nr. 




4 V 6 + 3 
which for 6=1 reduces to 



(b + 2) 
{N c + 2 



(38) 



(39) 



The results for q(NNir), o(AA7r) and g(ANn) have been 
obtained earlier 0, El Ej ■ Note the dependence of these 
quantities on N c + 2. In Ref. it was shown that quark 
matrix elements are given as an expansion in J 2 /(N c + 2), 



and the results Eq. I|35I36I38|I agree with the general form 
given there. 

The connection of the quark model to the soliton model 
is given by constructing spin-flavor coherent states for 
baryons containing N c quarks 



\V,N C 



(40) 



where the coefficient r\ is normalized to rf a rj* a = 1. The 
corresponding holomorphic wavefunction is 



(41) 



The chiral soliton in the standard configuration Uq, with 
(I + 3)Uq — corresponds to rf a = ?7q q , where 



(42) 



Rotated soliton configurations are given by applying the 
corresponding spin and flavor rotations to 77q°\ The inner 
product of two coherent states is 



(7?';iV c |ry;iV c ) = WO 



(43) 



In the limit iV c — > oo, this reduces to a (^-function in 77. 
The overlap of a soliton and a rotated soliton also van- 
ishes unless the rotation is the identity transformation, 
which is why the quark and soliton models give the same 
matrix elements as iV c — > 00 18]. The spin-flavor coher- 
ent states Eq. l|4T?|l were used in Ref. [T3 to compute the 
N c — ► 00 matrix elements in terms of simple integrals. 

Excited baryons such as the [70, 1~] are not in com- 
pletely symmetric spin-flavor state, since one of the 
quarks is in an orbitally excited state. The bosonic opera- 
tor formulation automatically constructs only completely 
symmetric states. One can include excited baryons by in- 
troducing a new operator q' for excited quarks. States 
generated by and q'^ have no symmetry restriction 
between the q^ and q't spin-flavor indices. This allows 
one to treat the [70, 1 _ ] baryons for arbitrary N c , using 
a core of N c — 1 ground state quarks, and one excited 
quark 0, H, 0, E3 • Similarly, one can treat heavy quark 
baryons by introducing operators c and b that act on 
c and 6 quarks, and constructing states of N c — 1 light 
quarks and one heavy quark |ll| . The advantage of using 
bosonic operators decreases rapidly with the number of 
quarks that must be treated as special. 
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